NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



DBG- 31426 


of Lagoos’s Orbit 

i!4SA-Ta-6073«) OM THE SECULAB DECREASE IN 
THE SEHIHAJOB AXIS OF LA6E0S OBBIT (NASA) 

43 p HC A03/HF AU1 CSCL 22A 

Duel as 

G3/13 32995 


David Parry Rufoincam 


JULY 1980 


National Aeronautics and 
Space Administration 

Goddard 8paco Flight Carter 

Greenbeit, Maryland 20771 



ON THE SECULAR DECREASE IN 
THE SEMIMAJOR AXIS OF LAGEOS’S ORBIT 

David Parry Rubincam 
Geodynamics Branch 
Code 921 


July 1980 


GODDARD SPACE FLIGHT CENTER 
Greenbelt, Maryland 20771 


CONTENTS 


Page 


ABSTRACT v 

1. INTRODUCTION 1 

2. RESONANCE 5 

3. GRAVITATIONAL RADIATION 7 

4. POYNTING-ROBERTSON EFFECT 8 

5. YARKOVSKY EFFECT 10 

6. SCHACH EFFECT 18 

7. TERRESTRIAL RADIATION 20 

8. MAGNETIC DESPIN OF LAGEOS 22 

9. DRAG FROM INTERPLANETARY DUST 23 

10. ATMOSPHERIC DRAG 24 

11. DISCUSSION 34 

ACKNOWLEDGMENTS 35 

REFERENCES 36 


?ing page bla:::; not 


FILMED 


iii 


ON THE SECULAR DECREASE IN 
THE SEMIMAJOR AXIS OF LAGEOS’S ORBIT 

David Parry Rubincam 
ABSTRACT 

The semimajor axis of the Lageos orbit is decreasing secularly at the rate of -1 . 1 mm day -1 due 
to an unknown force. Nine possible mechanisms are investigated here to discover which one, if any, 
might be the force. Five of the mechanisms, resonance with the earth’s gravitational field, gravita- 
tional radiation, the Poynting-Robertson effect, transfer of spin angular momentum to the orbital 
angular momentum, and drag from near-earth dust are ruled out because they are too small or require 
unacceptable assumptions to account for the observed rate. Three other mechanisms, the Yarkovsky 
effect, the Schach effect, and terrestrial radiation pressure could possibly give the proper order-of- 
magnitude for the decay rate, but the characteristic signatures of these perturbations do not agree 
with the observed secular decrease. Atmospheric drag from a combination of charged and neutral 
particles is the most likely cause for the orbital decay. This mechanism explains at least 71 percent 
of the observed rate of decrease of the semimajor axis. It probably explains all of the decay since 
(a) the estimate of charged particle drag is conservative and (b) there may be substantial quantities 
of neutral helium at Lageos ’s altitude, which helps to solve the “helium problem” accouting for its 
escape from the earth’s atmosphere. 
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ON THE SECULAR DECREASE IN 
THE SEMIMAJOR AXIS OF LAGEOS’S ORBIT 

1. INTRODUCTION 

Lageos, the Laser Geodynamics Satellite, was launched from the Western Test Range into a 
retrograde orbit on 4 May 1 976. It is a completely passive satellite with an outer shell of aluminum. 

Its spherical surface is studded with laser retroreflectors. Inside it has a cylindrical brass core. Lageos’s 
orbit is very nearly circular, with a semimajor axis of about 2R E . The orbital inclination with respect 
to the earth's equator is about 1 10 degrees. Other relevant data on Lageos are given in Table 1 and in 
Smith and Dunn (1980). 

Lageos’s orbit is suffering an acceleration which is not presently modeled in the orbit determina- 
tion computer program (Smith and Dunn, 1980). This acceleration is causing the semimajor axis to 
decrease secularly at the rate of -1 .27 X 1 0" 8 ms" 1 * -1 . 1 mm day" 1 (see Fig. 1 ). In an earlier short 
paper (Rubincam, 1980) we attributed this decay to atmospheric drag, either from the charged or 
neutral particles in the earth’s atmosphere or a combination of both. The object of the present paper 
is to discuss this conclusion in more detail, as well as show that other possible mechanisms which 
might explain the orbital decay appear to be inadequate to do so. These mechanisms are: resonance 
with the earth’s gravitational field, gravitational radiation, the Poynting-Robertson effect, the Yarkovsky 
effect, the Schach effect, terrestrial radiation, magnetic despin of the satellite, and drag from interplane- 
tary dust. The Schach effect is believed to be new and presented here for the first time. 

We will usually be content with order-of-magnitude arguments in our investigation. Also, we 
use SI units throughout, plus centimeters and millimeters when discussing lengths and days and years 
when discussing time. These latter non-SI units prove to be very convenient when considering the 
rate of decay of the semimajor axis. Most of the notation and numerical data are given in Tables 1 
and 2. Other symbols are explained as needed in the text. 
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Table 1 

Notation for quantitiei related to Lageot. The quantities in the first part of the table 
refer to Lageos’s orbit, while the second part refers to the satellite itself. Dash (-) in- 
dicates various numerical values are assumed. 


Quantity 


Symbol 


Numerical Value 


semimajor axis 


1.227 X 10 7 m 


eccentricity 


0.004 


true anomaly 


inclination 


109.9 deg 


mean anomaly 


mean motion 


4.65 X 10“* s 


-4 .-1 


velocity 
acceleration 
argument of perigee 
nodal position 


cross-sectional area 


0.2827 m 2 


drag coefficient 


411 kg 


radius 


0.30 m 


10.3 rad s -1 


temperature 
initial spin rate 











Table 2 

Notation for quantities related to the earth and sun (first part of the table) and to the uni- 
verse (the second part). Symbols not shown in Tables 1 and 2 are defined in the text. 

Dash (-) indicates various numerical values are assumed. 


Quantity 

Symbol 

Numerical Value 

gravity field coefficients 

^sm* 

- 

solar flux 

F. 

1.35 X 10 3 Jm“ 2s_1 

mass of earth 

M E 

5.97 X 10 24 kg 

radial distance from center of earth 

r 

- 

radius of earth 

r e 

6.371 X 10 8 m 

rotation angle for earth 

e 

- 

rotation speed of earth 

d 

7.29 x 10" 5 rad s" 1 

speed of light 

c 

3 x 10 8 ms" 1 

permittivity of space 

«o 

8.85 X 10 12 FnT 1 

gravitation constant 

G 

6.67 X 10" 11 m 3 kg -1 s" 2 

Boltzmann's constant 

k 

1.38 X 10" 23 J deg" 1 

proton charge 

q 

1.60 X 10" 19 C 

Stefan-Boltzmann constant 

0 

5.67 X 10" 8 W m" 2 deg" 4 

time 

t 

- 
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2. RESONANCE 


We flnt enquire into the possibility that the apparent secular decrease in the semimajor axis is 
due to a long-period resonance with the earth's gravitational field. 

Frm Kaula ( 1 966, p. 40) the change in the semimajor axis due to terms of degree 8 and order m 
in tite spherical harmonic expansion of the field is 

gm eY /2 / r eY 

a J \~Tj ^ +< ^®* m p<l 

<8-2p)w+(8-2p + q)& + m(6-0) 



where (Kaula, 1966, p. 37) 


smpq 


"Sm 


-S 


tin 


rs (m i 

L C BmJ 


e-m even 

s-m odd 
s-m even 

e-m odd 


cos ((8 - 2p)w + (8 - 2p + q) M + m (ft - 0)] 


sin 1(8 - 2p)w (8 - 2p + q) M + m (ft - 0) j . 


( 2 ) 


Here C 8m and S 8m are the unnormalized spherical harmonic coefficients and o>, M, 0, 0 are assumed 
to vary linearly with time. The denominator of ( 1 ) is nearly zero when the orbit is close to a resonance. 

We can simplify (1) in the following manner. Since (2) varies sinusoidally with time, it will have 
its maximum rate of change when we choose the origin of time t such that (2) becomes 


Ssmpq a C; m rin [(8-2p) W + (8-2p + q)M + m(fl-0)] 

* C' m sin {l(8 - 2p)w + (8 - 2p + q)M + m(6 - 0)j t) 

= C; m l(fi-2p)<j+(fi-2p+q)M + m(n-0)]t 

for small t, where C 8m is a typical coefficient of degree 8 and order m. Substitution of the above 
equation into (1) shows that the dotted expression in the numerator cancels with that in the de- 
nominator, leaving 


Aa 



W^tW^E^P+q) 



t 
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which varies linearly with time, to tint obviously the expression in curly brackets is 


to /GM e \i/ 2/R e V 

T, m2 {~) \ r] W> <W e) (« - Jp ♦ 1) 

where we have switched to normalized coefficients and inclination functions: 


( 3 ) 


'tm 


(£ + m)! 


1/2 


W *>' 


(2fi + l)(2-6 om )(fi-m)! 
■(2£+l)(2-5 om )(8-m)n‘/ 2 


C 

v «m 


^*mp^‘ 


(£ + m)! 

We now simplify (3) in order to estimate its magnitude. We first look at the eccentricity func* 


tionG gpq (e). Since 


G !p< j(e) 




e^ 1 , q W 0 
q-0 


the values for q not near zero are negligible due to the small eccentricity of the orbit. Hence we take 


q * 0 

£-2p + q*£-2p 
G ew (e)* 1. 

Next, we note that p runs from 0 to £, so that (4) becomes 

|£- 2pl < £. 


And since a « 2R E , we take 



Finally, from Kaula's (1966, p. 98) rule-of-thumb we use 


ic; 


em 1 


10“ 5 

«2 


(4) 

(5) 

(6) 

(7) 

( 8 ) 
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Substituting (4) - (8) in (3) and giving numerical values to GM E and a yield 


da 0-ii4F g mp(1) _j ■» ■ imp 


10? . -i 

m ,-» m mm day 1 . 


dt g . 2 « S • 2* 

All that remains is to estimate F gmp (l). Numerical experiments indicate that at most 


■W»l * »• 


Using this in the previous equation we have 

I da I 10 3 * * * 7 


. . » ' ■ mm day - * 

dti g. 2 * 


(9) 


as an estimate of the rate of change of the semimajor axis with time due to resonance with terms 
of degree fi. 

The rate given by (9) fails below the observed rate when £ * 22. Terms to degree and order 
36 have been numerically investigated to see if a resonance occurs with one of them to account 
for the observed rate. None of them do. Hence given (9) and the results of the numerical in* 
vestigation, it is unlikely that the observed rate of decrease of the semimajor axis is due to reso- 
nance with a term or even a number of terms above degree 36. We therefore conclude that reso- 
nance with the gravity field is probably not the unknown force operating on La geos. 


3. GRAVITATIONAL RADIATION 

The drag force easiest to rule out as the unknown force is gravitational radiation. Here the drag 
is due to radiation reaction; accelerated masses lose energy through gravity waves, just as accelerated 
charges lose energy through electromagnetic waves. The rate of change in the semimajor axis due to 
gravitational radiation is given by (Peters. 1964, cq. 5.6): 

d. _ 

dt *' Jc s,> 

where small terms involving the satellite mass and orbital eccentricity have been ignored. Substituting 
numerical values in the right side of the above equation yields 
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—• ft -1,4 x 10" 2 ® mm day -1 
dt 


(10) 


which fails to explain the observed decay by 26 orden of magnitude. Hence the perturbation is not 
due to gravitational radiation, which is of no importance for earth satellites. Orbital decay through 
this mechanism does, however, appear to have been observed in a binary pulsar system, providing an 
important test of general relativity (Taylor et al ., 1979). 

4. POYNTING-ROBERTSON EFFECT 

The Poynting- Robertson effect also produces a drag force (Robertson, 1937) It is known to be 
important for small particles in solar orbit (Stacey, 1977, p. 1 7-20; Lovell, 1954, pp. 402-409). It is 
accordingly investigated here for Lageos. 

The effect is due to the reradiation of light incident upon the satellite. Some of the light falling 
on Lageos is absorbed by it and is assumed to be reradiated isotropically in its own frame of reference. 
If the satellite is in motion with respect to an observer in another frame, such as the earth, then the 
observer sees a Doppler shift in the irradiated light. The light emitted in the direction of motion is 
shifted towards the blue end of the spectrum while the light emitted opposite to the direction of 
motion shows a red shit t . Since more energy and momentum is carried away from the satellite by 
the blue-shifted light than the red, the satellite feels a reaction force acting opposite to its direction 
of motion, it is this force which produces the drag. 

The light incident upon Lageos comes from two major sources: earth and sun. We investigate 
carthlight first. We assume that the earth radiates energy isotropically at the same rate at which it 
receives energy from the sun. The earth is also taken to be the stationary reference frame. 

Soter ctjd. (1977) gives a simple derivation for the Poynting-Robcrtson effect. The drag term 
in the force is 


m l v 


FA "* 

— V 


(II) 


where F is the integrated flux, giving a tangential disturbing function 
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in die Gaussian form of Lagrange’s planetary equations (Blanco and McCuskey, 1 961 , p, 1 78): 


da 

dt 



t(eninf)R + a(l -e 2 )SJ, 


(12 


to that 


da B 2S _ 2FA 
’ "nM L c 2 

We have assumed the orbit to be circular. Since F varies as r“ 2 . we take 

r Fq" 1 Ml 

r 2 a 2 

where F q is the flux at the earth's surface. Also, for circular ort.ts 


(13 


(14 


v- 

a l/2 


(15 


and by deftnition n is 


n ■ 


VGMg 

a 3/2 


(16 


Substituting ( 14) - ( 1 6) into ( 1 3) gives 

da ^ 2 FqR|A 
M t «*. ' 

We need to find F 0 in ( 1 7), the values of the other quantities being known. This is easily done 
the earth intercepts irR^F, power from the sun and reradiates 4»R|F 0 . Therefore 


where F, is taken from Robertson (1937) and shown in Table 2. Plugging numerical values in (17) 
gives 

-1.69 X 10" 11 ms" 1 
dt 

* - 1 cm per 1 8 yr (18) 

The Poy ntin g-Robertson effect due to earthlight is much too small to account for the observed rate 
of decrease in a. 

Sunlight is also too small to explain the observed orbital decay. Repeating the calculation using 
the direct solar flux F, in (13) gives a rate of decay 16 times larger than that due to earthlight (a fac- 
tor of 4 coming from F # = 4F 0 and another factor of 4 due to the r" 2 behavior of earthlight given by 
(14)), or 

— a -2.70 X 10" 10 ms" 1 
dt 


= - 1 mm per 43 days (19) 

The calculation is only approximate since (1 1) is valid for an inertial frame in which the source of 
light is at rest. The earth orbiting about the sun is neither an inertial frame, nor is the sun at rest in 
it. However, these corrections are small (Allan, 1967, pp. 74-75), as are the variations in flux due to 
the changing satellite-sun distance. The most serious neglect is the effect of the earth’s shadow. But 
the shadow merely turns the drag force off over part of the orbit, slowing the overall decay rate, so 
that it makes the Poynting-Robertson effect smaller than given in (19). We must conclude that the 
Poynting- Robertson effect is not the unknown force; it is too small. 


5. YARKOVSKY EFFECT 

We consider now the Yarkovsky effect (6pik, 1951, pp. 194-197; Lovell, 1954, pp. 41041 1). 
This is a differential thermal effect caused by light falling on a rotating satellite with appreciable 
thermal inertia. 
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Radiation from a light source warms the surface of the satellite. The hottest part of the satellite 
would be the sub-source point if the satellite were not rotating. However, suppose light is falling on 
a rotating satellite as shown in Fig. 2. Here areas on the satellite's surface are carried around from 
the shadowed side into the light and are wanned by it. There is a delay in heating up due to the 
thermal inertia of the satellite, so that the hottest part is the “afternoon" side of the satellite and not 
the sub-source point. This is similar to temperature variations here on earth, were afternoon is the 
hottest time of day, instead of noon. 

The asymmetric heating causes a net force to act on the satellite, since photons from the hotter 
areas of the surface carry away more momentum than the colder areas. This force, which makes a 
definite angle with the satellite-source line, perturbs the orbit. 

Light for the Yarkovsky effect comes from the earth and sun, just as for the Poynting-Robertson 
effect. We again investigate earthlight first. 

To obtain a qualitative idea of how the Yarkovsky effect from earthlight can affect the orbit 
we refer to Fig. 3. Here the geometry of the orbit and Lageos spin axis is such that there is a com- 
ponent of force acting consistently opposite to the satellite velocity vector. This component will of 
course cause the orbit to decay. If the satellite is spinning in the sense opposite to that shown in 
Fig. 3, then the force will increase the size of the orbit. Thus the effect on the orbit depends crucially 
on the orientation of the spin axis in space. The presumption at this stage of the investigation must 
be that the spin axis is oriented largely as shown in Fig. 3, in order to explain the observed decay 
shown in Fig. 1 . 

We tum now to estimating the magnitude of the acceleration due to the Yarkovsky effect. We 
will assume that Lageos is composed of two hemispheres, one at temperature T and the other at tem- 
perature T + AT. If each hemisphere radiates like a blackbody over a characteristic area of irR^, then 
the acceleration is 



oirRj[(T + AT) 4 -T 4 ] 


4*R{oT 3 AT 


( 20 ) 
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V 


ACCELERATION 

Figure 2. Schematic diagram illustrating the Yarkovsky effect on a satellite. Light im- 
pinges on the satellite from the top, as illustrated by the parallel vertical arrows. Some of 
the light is absorbed by the satellite and re-emitted. The re-emitted light is denoted by the 
arrows pointing outwanis from the satellite. The radiation reaction forces from this light 
give a net acceleration V to the satellite, shown as the thick arrow. The acceleration is not 
in the direction of the incident light, due to the rotation of the satellite. 
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assuming AT « T. Hie temperature T of a blackbody at distance a from the earth is (Bianco and 
McCuskey, 1961, p. 53): 




where T E a 280 K is the blackbody temperature of the earth. Substituting numerical values in (20) 
gives 

| V| = 1.43 X 10" 12 AT ms" 2 . (21) 


Taking 



we find 


Ud * 

— = -0.51 AT mm day -1 . (22) 

dt 

A temperature difference of 2 K in our model will explain the entire perturbation. Since this 
does not seem unreasonable, it would appear a detailed treatment of the problem involving the 
albedo of the satellite, its thermal properties, etc. is needed before we can make a judgment as to the 
importance of the Yarkovsky effect. However, we will leave our result (22) in abeyance for the mo- 
ment and show later that such a detailed treatment is not needed. We turn our attention now to the 
Yarkovsky effect due to sunlight. 

We can compute the approximate magnitude of the acceleration from (20), this time taking 
T = T E , since Lageos is at about the same distance from the sun as it is the earth. This yields in 
analogy to (21) 

I V I = 1.14 X 10 -11 AT ms' 2 (23) 

— ► 

To compute the effect of this acceleration on the orbit, let V have the components V = (f x ,f y ,f z ) 
in the coordinate system shown in Fig. 4. The orthogonal accelerations R, S, and W are then given by 


14 


R 


cos(a> + f)cosf2 

cos(w + f) sinft 

sin(w + f)sinl 





-cosIsinfisin(u + 0 

+coslcosf2sin(w+ f) 




S 

■ 

-sin(iu 4- f)cos SI 
-cosIsinftcos(w + 0 

-sin(cj+ f)sinf2 
+coslcosf2cos(w + f) 

cosfcc + f)sinl 


f y 

W 


sinking 

-sinlcosfl 

cost 


A 


where the rotation matrix comes from Goldstein (1950, p. 109). Substituting the expressions R and 
S into (12), noting that 

a(l-e 2 ) 

1 + e cos f 

and integrating the true anomaly f from 0 to 2ir yields the average rate of change for the semimajor 
axis 


da 

dt 


2ea 3 / 2 

y/m^ 


[(-sincocosfl - cos!sinJ2cosco)f x 


+ (-sinwsinfi + cosIcosf2cosoj)f y 

+(coswsinI)fJ (24) 

This procedure assumes that the magnitude and direction of V does not change appreciably over one 
revolution (certainly a good assumption), and that the orbit is oriented so that the satellite does not 
enter the earth’s shadow (to be discussed shortly). 

From (23) and (24) we find that the rate of change of semimajor axis with time is roughly 


da 

dt 


a 1.96 X 10 -10 AT ms -1 = 0.017 AT mm day -1 . 


(25) 


An unreasonably large temperature difference of 10 K is stilt much too small to contribute appreciably 
to the rate of change of the semimajor axis. Thus the Yarkovsky effect when the orbit is in full sun- 
light is unimportant. 
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7 

Figure 4. The orbital geometry used in computing the 
Yarkovsky effect from sunlight. 
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This remit does not hold if part of the orbit falls in the earth’s shadow. In this case V is zero in 
the shadow, which is equivalent to integrating f over less than lit when averaging the perturbation 
over one revolution. In fact, as much as one-sixth of the orbit call fall in the earth’s shadow. In these 
circumstances large terms in the expression for da/dt do not average out to zero and we have 
approximately 


da 

dt 


2a 3/2 ^ 

■- ■ I V I a 4.91 X 10" 8 AT ms' 1 

\/6m f 


= 4.24 AT mm day" 1 


(26) 


or a factor of e" 1 larger than given by (25) when the orbit is in full sunlight. This is even larger than 
(22), assuming the same temperature difference. We cannot rule out the Yarkovsky effect from sun- 
light using order-of-magnitude arguments only when the orbit is shadowed, just as we could not rule 
out the Yarkovsky effect from earthlight. 

We will, however, attempt to rule out the Yarkovsky effect from both sources on the basis of 
the orbital geometry and the known behavior of the semimajor axis shown in Fig. 1. 

Apogee kick stage separation occurred at +4.5 deg geocentric latitude and +20.8 deg east longi- 
tude. At this time the Lageos spin axis had elevation -1 1 .0 deg and azimuth +1 58.7 deg. The nodal 
position was +28.7 deg. The use of these data plus the known orbital inclination show that the spin 
axis was approximately in the plane of the orbit, as might be expected. The spin axis must have 
slowly moved out of the orbital plane as the node progressed along the equator, assuming the posi- 
tion of the spin axis stayed fixed in space. 

The Yarkovsky effect from sunlight or earthlight will have no long-term effect on the semimajor 
axis when the spin axis lies in the orbital plane, as may be seen from considerations like those shown 
in Fig. 3. Hence, in the beginning part of the curve shown in Fig. 1 , when the satellite is in full sun- 
light (so that the contribution from the sun is negligible) and the spin axis is nearly in the orbital 
plane (so that the contribution from the earth is negligible), the curve should have near zero slope, 
assuming the unknown force is due to the Yarkovsky effect. The slope should become steeper as 
the spin axis moves out of the orbital plane and the contribution from the earth becomes appreciable. 
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This don not occur, u is obvious from Fig. 1 . If anything the dope is steeper at this part of the 
curve than at later times and then flattens out a little. This is just the opposite from what we would 
expect from the Yarkovsky effect. We conclude that the unknown force is not due to the Yarkovsky 
effect. Moreover, the temperature differences across Lageos must be small, since the curve deviates 
but little from a straight line, indicating the Yarkovsky effect from both the earth and sun are operat- 
ing at a low level. Still, analysis of the irregularities in the curve may yield information about the 
Yarkovsky effect. 

We should note that we have assumed two different temperatures for Lageos: one due to sun- 
light and the other to earthlight. This inconsistency makes no essential change in our argument, since 
we were interested only in obtaining a rough estimate of the possible importance of the Yarkovsky 
effect. 

6. SCHACH EFFECT 

Related to the Yarkovsky effect is the Schach effect, which was discovered by Milton Schach of 
the NASA Goddard Space Right Center in the course of this investigation. Its operation depends 
upon the coding off and heating up of the satellite as it moves in and out of the earth’s shadow. We 
refer to Fig. 5 for its explanation. 

The satellite is wanned up by sunlight when it is outside of the shadow. When it enters the 
shadow it will take some time to cool off, due to its thermal inertia. Likewise, it will take some 
time to warm up when the satellite exits the shadow into sunlight. Qualitatively we can think of the 
satellite cooling off and heating up instantaneously as it moves through a rotated shadow, as indicated 
in Fig. 5. 

This can affect the orbit in two ways. One is a delay in turning off and on the Poynting- 
Robertson effect due to sunlight. This causes no essential change in the long-term operation of the 
Poynting-Robertson effect or our arguments based upon it. The other way it can affect the orbit is 
when the satellite is a poor heat conductor and rotating slowly (or not at all) or with its spin axis 
pointing towards the sun. In this case the hemisphere facing the sun becomes hot, so that the 
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Figure 5. Schematic illustration of the Schach effect. A satellite is shown at two points 
along its orbit. Note that the acceleration $ is in the direction of the sunlight (so that the 
Yarkovsky rotation is ignored) when the satellite is not in shadow. If the satellite cools off 
and warms up instantly when moving through the shadow, then the acceleration is zero in 
die shadow and the acceleration over arc FED tends to cancel with that over arc FAB by 
symmetry. However, if the satellite has some thermal inertia, then it takes some time to 
cod off and heat up when moving through the shadow. Qualitatively we can think of a 
thermally inertialess satellite passing through a rotated shadow inside of which t? is zero. 
In this case the acceleration over arc FE cancels that over arc FA, but no cancellation occurs 
for arc ABC. This leads to a net acceleration when averaging over one revolution which 
tends to increase the semimajor axis. 
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acceleration due to the reradiated sunlight points in the direction opposite to the sun (ignoring the 
Yarkovsky effect). This acceleration does not cancel over all parts of the orbit, giving a net impulse 
to the satellite when it is averaged over one revolution, as explained in the figure caption. 

The Schach effect cannot be the unknown force causing the orbit to decay. For one thing, the 
force operates whether the orbit intersects the shadow or not, while the Schach effect operates only 
when La geos moves through the shadow. For another the force acts to increase the semimajor axis, 
as should be clear from Fig. S. Further, the Schach effect must operate at a low level, like the 
Yarkovsky effect, since the curve shows only small deviations from a straight line. Again the analy- 
sis of the deviations may be of some interest in trying to detect the effect. 

While the Schach effect is not important for Lageos, it may have some interesting consequences 
for the orbits of the particles comprising the rings of Jupiter, Saturn, Uranus, and possible microtek- 
tite rings of the earth (O’Keefe. 1980). But we will not pursue this here, being off the main topic. 

7. TERRESTRIAL RADIATION 

Yet another radiation force arises fromthe pressure of terrestrial radiation on the satellite. This 
topic has been treated by many authors, including Wyatt (1963), Sehnal (1970), Prior (1970). and 
Smith (1970). Recent work Ms been done by Saslaw (1978. 1979) and Lautman (1977a, 1977b). 

See the references of these papers for other relevant work. 

Terrestrial radiation falling on the satellite comes from three sources (Smith. 1 970). One is the 
diffuse reflection of sunlight from land, water, clouds, etc. Another is the specular reflection of sun- 
light from still water, ice, or snow. These two together constitute the “albedo problem,” although 
most attention is given to diffuse reflection. Last is the thermal infrared radiation coming from the 
absorption of sunlight and its reradiation by the earth. 

Analytical treatment of diffuse reflection is difficult because of the complicated way it varies 
in both space and time. The simple early models of diffuse reflection are probably unrealistic. Hence, 
to obtain a general indication of how diffuse radiation affects a satellite, we will rely on the work of 
Lautman (1977a, 1977b), who gives the most detailed treatment of the problem to date. 
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Lautman (1977b) assume* the albedo a, varies according to the equation a s ■ a 0 + a 2 sin 2 d, 
where $ to latitude and a 0 , a 2 are constants. It to also assumed that “the perturbing acceleration is 
given by the nonterminator expressions during half of the mbit that the satellite to closest to the sun 
and that it to zero during the other half,” with the more complicated terminator expressions to be 
given in a future paper (Lautman, 1977b, p. 3). He finds (see his Fig. 5) in a sample calculation for 
the Page os balloon satellite that the semimajor axis a decreases by about 160 m over 24 days and 
then increases by about 280 m over the next 36 days and was still increasing when the 60 day inte- 
gration time was reached. The following parameters were assumed in this calculation: a * 1.06 x IQ 7 m. 
e- 0.0628, 1 - 86.9 deg, ft * 329.8 deg, w • 215.8 deg, T e - MJD 39384.0, a 0 - 0.219, a 2 - 0.410, 
and the area-to-mass ratio was 1 1 .8 m 2 kg -1 . 

Both Pageos and Lageos are in polar orbits; their semimajor axes are comparable, and the eccen- 
tricities of their orbits are small. Hence it is reasonable to assume that the qualitative behavior of 
Lageos’s orbit will be similar to Pageos *s. This qualitative behavior is inconsistent with what is ob- 
served on Lageos. Pageos’s semimajor axis decreases and then increases over a short time period of 
60 days. Lageos’s semimajor axis shows a steady decrease over the long time period of 3 years. Since 
diffuse radiation from the earth is capable of increasing the semimajor axis over short time intervals 
(Hit Lageos’s semimajor axis docs not do so, we conclude that diffuse radiation is a most unlikely 
candidate for being the unknown force acting on Lageos. 

It is of some interest to estimate the inaptitude of the perturbation due to diffuse reflection 
for Lageos. If Lageos were in Pageos's orbit, then the change in a would be very roughly 1.6 cm over 
the time involved. This is due to the very much smaller area-to-mass ratio for Lageos. If the orbit 
had Lageos's eccentricity of 0.004, then the change in a would be about 16 times smaller, or about 
1 mm over the time involved, since the long-period terms in Lautman’s (1977b) expression for the 
change in a (his equation 34) are proportional to the first or greater power of e. Thus it would ap- 
pear that the magnitude of the diffuse reflection effect is small compared to that of the unknown 
force is small. This to consistent with the small magnitude of the fluctuations about a straight line 
shown in Fig. 1 . 


21 


Wyatt (1963) discussed the specular component of albedo radiation pressure. Again the mathe- 
matical difficulties in producing a realistic model are great, but a rough estimate of its magnitude can 
be made. Wyatt*! (1963) work suggests that the specular component is probably only a few percent 
the size of the diffuse component. The latter is already judged to be small. Further, we would ex- 
pect its operation to be erratic, since the necessary conditions of calm water or snow or ice fields and 
cloudless skies above them are erratic. This is inconsistent with the regularity of behavior shown in 
Fig. 1. We dismiss the possibility of specular reflection being the unknown perturbation for these 
reasons. 

The effect of infrared radiation pressure on Lageos is probably also too small. We would expect 
the intensity of the radiation to be roughly that of diffuse reflection, but more symmetrically dis- 
posed over the earth. Since the effect of diffuse radiation is judged to be small, and since symmetry 
tends to produce no secular effects as found by Wyatt (1963), we will not consider it further. 


8. MAGNETIC DESPIN OF LAGEOS 

Lageos 's spin rate is being slowed down by the earth's magnetic field via eddy currents (Zonov, 
1961 ; Smythe, 1950, pp. 390420). If we assume that the lost sp.n angular momentum is gained by 
the orbital angular momentum in analogy to tidal friction, then this will affect the semimajor axis. 
We will compute the magnitude of this effect. 

Assume a circular orbit; the orbital angular momentum L is then (Blanco and McCuskey, 1961, 
p. 133): 

L- % /ES^M L a W 


so that the change in the semimajor axis Aa is then 

2» ,/2 AL 


Aa 1 


v^GMg M l 


(27) 


due to a change in angular momentum AL. assuming the spin axis is perpendicular to the orbital 
plane. We take AL to be the entire spin angular momentum of the satellite 
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|AL| -Cw L 


where we hike the moment of inertia C to be that of a homogeneoui sphere: 

C-0,4M l R[. 

Using the two equations above in (27) and numerical values from Tables I and 2 yield 

| Aa| w 1.28 x I O’ 4 m« 0.128 mm. 

Since this is the entire effect of stopping the spin and not just a rate, there is no possibility of mag- 
netic despin being the unknown force acting on Lageos's orbit. 


9. DRAG FROM INTERPLANETARY DUST 


Collisions with dust particles will cause a drag force to be exerted on Lageos. We will use the 


standard drag equation (Blanco and McCuskey, 1961, p. 204) 

j. _ c„A,vi 

dt M L n 


(28) 


to estimate the magnitude of the drag force, where in this equation p is the density of near earth dust 
in kg m~ 3 and a circular orbit has been assumed. Using ( 1 5) and (16) in (28) and numerical values 
from Tables 1 and 2 give 


da 

dt 


CpAp^ONgi 


* -4.81 x !0 7 Corns' 1 


(29) 


for Lageos. We now need C D and p. 

For C D we take the standard value 2.2 (Cook, 1965, p. 929). From Hughes ( 1975) we discover 
that the upper limit on the interplanetary dust density p ( at 1 AU from the sun is 

p,» 1.5 x 10- ,9 kgm- 3 . 

The absence of dust belts around the earth indicates that a generous estimate in the concentration of 
dust near the earth is (Shapiro, etal., 1966): 

pm I0pj« 1.5 X 10* 11 kgm" 3 
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Using these values for C D and p in (29) yield 



—a -1.59 X KT 10 ms -1 = -0.0137 mm day -1 (30) 

for the rate of decrease in a due to dust drag. This is two orders-of-magnitude too small to explain 
the observed rate. Since our assumptions were rather generous in making this estimate, we conclude 
there is no possibility of attributing the orbital decay to this mechanism. 

10. ATMOSPHERIC DRAG 

We finally consider atmospheric drag from neutral and charged particles as being the unknown 
force acting on Lageos. We look at neutral particle drag first. We can rewrite (29) as 

da , . , 

mC D N(a) = 1.58 X 10 11 m" 3 (31) 

Av/GM^a dt 

where the right side of (30) uses the observed decay rate of -1.27 X 10“ 8 ms -1 , and where p=N(a)m, 
with N(a) being the number density in particles nC" at Lageos’s altitude and m being the particle 
mass. This equation must be satisfied if we are to ascribe the observed orbital decay as being solely 
due to some constituent of the neutral atmosphere. 

We confine the discussion to hydrogen and helium only. The atomic hydrogen comes from the 
photodissociation of water vapor, while helium comes from the decay of uranium and thorium inside 
the earth. Other constituents of the earth's atmosphere, such as nitrogen, are too heavy to be found 
in appreciable amounts at Lageos 's altitude (Jacchia, 1977, Table 10). 

We consider hydrogen first. We take m in (3 1 ) as the mass of the hydrogen atom. Cook ( 1 965, 
Table 2) finds that C D tends to increase slowly with height, ranging from 3.6 at 1 000 km altitude to 
3.8 at 3000 km altitude for an exospheric temperature of 1 000 K. We will assume C D = 3.8 for 
hydrogen at Lageos’s altitude for all temperatures. In this case (31 ) reduces to 

N(a) = 4.15 x 10 10 m -3 (32) 

which is the number density necessary to explain the observed orbital decay due solely to neutral 
hydrogen drag. We now need to estimate N(a) for hydrogen and compare it to (32). 
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We estimate N(a) using the theory of planetary coronae (Chamberlain 1963, 1968; see also 
Brinkmann, 1971 and Chamberlain and Smith, 1971). For a spherically symmetric, nonrotating 
atmosphere the number density at radial distance r from the center of the earth is given by Chamber- 
lain (1963, p. 906): 

N(r) = N(r c )e" (Xc “ X) f(X) (33) 

where 

GM E m 

X = X(r) = — — . (34) 

kT c r 

Here r c and T c are the radial distance to, and temperature at, the critical level, respectively, while 
X c = X(r c ). The critical level is the altitude at which a neutral particle escaping towards infinity has 
something like a 1 /e chance of colliding with another molecule, where e is the base of the natural 
logarithm (Blanco and McCuskey, 1961 , p. 58). We take the critical level to be 600 km altitude, so 
that r c = 6371 + 600 = 6971 km. The f(X) appearing in (33) is a partition function which in turn is 
the sum of three other functions: f(X) = f bal + Numerical values for the f bal , and f esc 

for various values of X may be found in Chamberlain’s (1963) Tables 1 and 2. The values for N(r c ) 
at temperature T c come from Jacchia’s (1977) Table 10. 

Table 3 shows the number density for atomic hydrogen for various temperatures T c for Lageos’s 
altitude using (33) and (34). It also compares the number densities at 2500 km altitude using (33) 
and (34) with those of Jacchia’s (1977) model. The agreement between the two is quite good. 

Table 3 shows that N(a) reaches a maximum of about 0.54 X 1 0 10 m -3 near T c = 1 000 K at 
Lageos’s altitude. Thus the number density found from the theory of planetary coronae falls short 
by at least a factor of 8 from the number required by (32) to explain the orbital decay as being due 
entirely to neutral hydrogen, regardless of exospheric temperature. Hence unless the number density 
at the critical level or the drag coefficient are greatly increased somehow, it appears that neutral 
hydrogen drag can explain no more than about 12 percent of the observed decrease in a. 

We turn now to drag from neutral helium. Only 4 He need be considered, since the abundance 
3 He is too small to exert appreciable drag on Lageos (Mac Donald, 1963). We find 
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Table 3 

Number densities for hydrogen at various exospheric temperatures. The first column 
shows the temperature. The second shows the number density at 2500 km altitude 
computed according to equations (33) and (34), and the third column the number 
density at 2500 km according to Jacchia’s (1977) Table 10. Note that the agreement 
between the two is quite good. Hie last column shows the number density at Lageos’s 
altitude. 


T C (K) 

N(10 10 m" 3 

2500 km 
(computed) 

2500 km 
(Jacchia) 

5856 km 

(computed, Lageos) 

600 

4.28 

4.30 

0.35 

800 

3.26 

3.27 

0.50 

1000 

2.44 

2.45 

0.54 

1200 

1.88 

1.87 

0.54 

1400 

1.48 

1.47 

0.50 

1600 

1.19 

1.17 

0.46 

1800 

0.97 

0.96 

0.42 

2000 

0.81 

0.79 

0.38 

2200 

0.69 

0.67 

0.35 

2400 

0.59 

0.57 

0.31 

2600 

0.51 

0.49 

0.29 
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- b jnfc 









C D N(a) ■ 3.95 X 10 10 m -3 

from (31) where the mass m appearing in that equation is now 4 times that of hydrogen. 

Cook (1965, Table 2) finds file values for C D for helium to be somewhat lower than for hydro- 
gen, ranging from 2.9 at 2000 km altitude to 3.0 at 3000 km altitude for a temperature of 1 500 K. 
He also finds that C D * 3.1 for the same two altitudes for a temperature of 2000 K, Since C D is 
rather insensitive to both altitude and temperature, we will set C D * 3.0 without any great error. 

The above equation now becomes 

N(a) = 1.32 X 10 10 m -3 . (35) 

This is the number density necessary to explain the secular decrease in a as being due solely to neu- 
tral helium drag. 

The exospheric temperature needed to produce the number density (35) is 2200 K. The varia- 
tion in exospheric temperature since Lageos’s launch is shown in Figure 6. The curve is based on the 
solar flux and geomagnetic activity in the following manner (Jacchia, 1977, equations 20, 31a, and 
31b): 

T c = 5.48 F 0 - 8 + 101.8 F 04 + A g T c 

where 

* g T c = / A p si " 4 d0 i 

and 

Ap = 57.5 K p [ 1 + 0.027 exp (0.4 K p )] . 

Here F is the solar flux in 10 4 Jansky (10 -22 W m -2 Hz -1 ), F is t'ie average of F over six rotations, 
is the invariant magnetic latitude, and K p is the geomagnetic activity index. Obviously the tem- 
perature at no time approaches the 2200 K necessary to explain the orbital decay as helium drag. In 
fact, assuming a typical temperature of 1000 K from Fig. 6 in (33) and (34) gives N(a) = 10 7 m" 3 . 
The number density thus derived from the actual exospheric temperature falls about 3 orders-of- 
magnitude short of the necessary number (35). Neutral helium appears to be a poor candidate for 
explaining the secular decrease in a. 
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Figure 6. The exospheric temperatures as a function of time since Lageos’s launch, 
based on the solar flux and geomagnetic activity. 



In spite of this, however, we do not rule it out, since helium drag helps solve the “helium 
problem.” 

The helium problem is reviewed by MacDonald (1963), Kockarts (1973), Hunten (1973), Hartle 
and Mayr (1976), and Chamberlain (1978). It consists of accounting for the outgassing from the 
earth’s interior, escape, and abundance of both 3 He and 4 He in the earth’s atmosphere. The influx 
of 4 He from radioactive decay greatly exceeds the rate at which it is being lost through known 
mechanisms. A loss mechanism must exist, since otherwise its present abundance would have been 
built up in only 10 6 years. Further, if one assumes the loss mechanism for 4 He is thermal escape 
(as considered here), with the temperature of the exosphere being elevated in some unknown manner, 
then the loss of 3 He cannot be accounted for by thermal escape (Chamberlain, 1978, p. 278). As 
Chamberlain (1978, p. 278) says, “The problem will not go away and is unsolved.” 

The exospheric temperature required for the thermal escape flux of 4 He to be equal to its in- 
flux into the atmosphere from the earth’s interior was found by MacDonald (1963) to be 2200 K. 
This is the same temperature required to explain the drag on Lageos as being due entirely to helium. 
Hence helium is an attractive candidate for being the drag on Lageos, since it takes care of the drag 
force and the 4 He part of the helium problem in one stroke. However, if this mechanism is at work 
then the temperature of the exosphere would have to be elevated by an as yet unrecognized heat 
source, as pointed out by MacDonald (1963). Moreover, it would have to be a steady source largely 
uncorrelated with the solar flux and geomagnetic activity, as demanded by the lack of correlation 
between Figs. 1 and 6. Finally, it does not help the 3 He part of the helium problem. 

We wish to point out that our assumption of a spherically symmetric atmosphere is oversimpli- 
fied. The helium distribution shows a distinct bulge over the winter hemisphere (Kockarts, 1973) 
which has recently been directly measured (Mauersberger et al ., 1976). 

The thermal mechanism is not the only one which can increase the helium abundance in the 
upper atmosphere. Winds can greatly enhance helium escape, as shown by Hartle and Mayr (1976). 
Also, helium may become ionized and escape over the poles on the open magnetic field lines 


(Chamberlain, 1978, p. 278). There is the further possibility that charge exchange between He or H 
and hot He 4 in the plasm asphere produces energetic He atoms (Chamberlain, 1978, p. 278). It is of 
some interest to note that die trend of both theory and measurement has been towards a higher 
abundance of helium in the upper atmosphere (Nicolet, 1961; Young et al ., 1 977 ; Hartle and Mayr, 
1976). And, as pointed out before, an escape mechanism must exist; a high helium density at 
Lageos:$ altitude may allow a steady-state between influx and loss, removing the need for assuming 
that we live in a special period in which the helium abundance is building up (Kockarts, 1973). 

Helium may well play a significant role in the drag on Lageos. 

We turn our attention now to charged particle drag. A satellite can become charged from col- 
lisions with charged particles and from the photoejection of electrons from the satellite’s surface. 

The satellite charge will interact electromagnetically with the charged particles in its vicinity and 
cause the satellite to lose momentum; whence the drag (Chopra, 1961 ). 

Let us first estimate the number density of the charged particles. The dominant ions at Lageos’s 
altitude are H 4 , He + , and 0 + with number densities of about 2 X 10 9 , 2 X 10 8 , and 2 X 10 8 m” 3 , 
respectively. (Chappell et al ., 1970, p. 51). In Rubincam (1980) we assumed that H 4 was the only 
ion present with a number density of 10* m -3 . Thus the data of Chappell etah (1970) show that 
the charged particle density at Lageos’s altitude is much greater than assumed in the previous paper. 

A conservative estimate of the drag due to H 4 , He 4 , and 0 4 can be made by assuming that these 
ions behave like neutral atoms. In other words, we take the drag coefficient C D to be 3.8 for H 4 and 
3.0 for He 4 , just as we did in computing the neutral particle drag for these constituents. We take C D 
to be equal to the standard value of 2.2 for 0 4 (Cook, 1965, p. 943). 

The contribution of each type of ion to the decrease in the semimajor axis using the above values 
for number densities and drag coefficients can be computed from (29). The results are summarized in 
Table 4. H 4 can account for about 4.4 percent of the observed decay rate, while He 4 and 0 4 account 
for 1.4 and 4.1 percent, respectively. Altogether the ions account for about 10 percent of the ob- 
served decay rate, assuming the ions drag like neutral atoms. The neutral atmosphere accounts for 
another 10 percent (Table 3) as explained previously, leaving aside the problematical contribution 
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Table 4 


Amount of observed drag explained by charged particle drag, assuming the ions 
behave like neutral particles. The drag due to neutral hydrogen is also given. 


Constituent 

Number Density N(a) 
(m-3) 

Drag CoefTicient 
C D 

Percent of 
Observed Drag 

H+ 

2.0 X 10 9 

3.8 

4.4 

He + 

2.0 X 10® 

3.0 

1.4 

0 + 

2.0 X 10 s 

2.2 

4.1 




Ion subtotal: 9.9 

Neutral H 

4.5 X 10 9 

3.8 

Total: 10.0 




19.9 










of neutral helium. Thus we can conservatively attribute at least 20 percent of the observed rate of 
decrease in the semimajor axis of Lageos to atmospheric drag from the combined effects of charged 
and neutral particles. 

Charged particle drag is in general greater than neutral particle drag for equal particle densities 
(Chopra, 1961). Unfortunately, the various theories developed to assess the importance of charged 
particle drag on satellites do not agree on how much greater. For example, Chopra ( 1 96 1 ) in his 
Table VI compares values for C D found from the Jastrow-Pearse, Kraus-Watson, and Chopra-Singer 
theories of Coulomb drag. The values range from a low of about 2 in the Kraus-Watson theory, to 
about 6 in the Jastrow-Pearse theory, to a high of about 8 X 10 3 in the Chopra-Singer theory for a 
high altitude (20,000 km) satellite charged at -21 v. In Rubincam (1980) we used the Chopra-Singer 
result to show that the decay of Lageos’s orbit falls within the limit for drag predicted by this theory. 

However, the experiments by Knechtel and Pitts (1964) using mercury ions to measure charged 
particle drag in the laboratory indicates that the Jastrow-Pearse theory is more nearly correct. So we 
will use their results to estimate the importance of charged particle drag on Lageos. 

Knechtel and Pitts (1964) find that multiplying the charged particle drag for an uncharged satel- 
lite by 

q0 

s — 1 - 1.33 ~ (36) 

E 


gives the total (neutral plus charged particle) drag on the satellite for R L /X D ® 7.5, where q is the 
proton charge, E is the kinetic energy of the ions, <t> 0 is the satellite potential, and \ D is the Debye 
length (Reitz and Milford, 1967, p. 272): 


,kT 


Xn " 


D V 


2N(a)q : 


T 


(37) 


where T is the temperature of the ions. Equation (36) is derived from Fig. 5 of Knechtel and Pitts 
(1964). 

Since Lageos orbits inside the plasmasphere where the thermal regime is one of cold ions (Rus- 
sell, 1972), we will take T • 2000 K and E * 3kT/2. The number densities N(a) for H + , He + , and 0 + 


have already been given. Using these data in (37) give R L /X D = 6.15 for H + , and 1 .94 for He + and 
0 + . These values are lower than the R L fk D * 7.5 required by (36), but Knechtel and Pitts ( 1 964) 
did not run the experiment for these low values. However, since the factor s increases as R L /X D de- 
creases (Knechtel and Pitts, 1964, Fig. 5), we can use (36) to estimate a lower bound on the contribu- 
tion of charged particles to the drag. In this case (36) becomes 

s = 1 - 5.14 0 O (38) 

It remains to And the satellite potential 0 O to use in (38). Satellites at the geosynchronous 
altitude can become charged up to hundreds or even thousands of volts (e.g. Garrett et al ., ( 1 980). 

At lower altitudes (~ 1000 km) satellite potentials appear to be on the order of a volt (Samir et al. . 
1979, p. 102). We will assume the lower figure is more nearly correct and take <f> 0 = - 1 v. (The 
voltage is negative si nee electrons collide with a satellite much more frequently than with positive 
ions, because of their higher speeds at any given temperature.) 

Using this voltage in (38) yields s * 6.14. We will approximate the ion drag for an uncharged 
satellite by the neutral particle assumption (Table 4)). We have already found that the ions account 
for about 10 percent of the observed decrease in the semimajor axis if they behave like neutral par- 
ticles. Thus multiplication of this figure by s gives 61 percent as the amount of the observed decay 
which can be explained as being due to charged particle drag, using the data of Knechtel and Pitts 
(1964). Neutral hydrogen drag accounts for another 10 percent, so that 71 percent of the observed 
drag can be reasonably attributed to the combined effects of charged and neutral particle atmospheric 
drag. A similar conclusion was reached by Barlier (1980). 

We conclude our investigation of atmospheric drag by saying that if we can reasonably attribute 
71 percent of the observed decay to a combination of neutral and charged particle drag, then it is 
likely that this mechanism explains the entire effect. This is because our estimate of the magnitude 
of charged particle drag is a probable lower bound and because the 71 percent figure does not include 
a possible neutral helium contribution. 
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11. DISCUSSION 


We have examined nine mechaniams to discover which one, if any, might be die force causing 
the semimgjor axis to decrease at the rate of -1 .1 mm day’ 1 . Gravitational radiation and transfer of 
spin angular momentum to orbital angular momentum are incontrovertibly too small to be the un- 
known force. The Poynting- Robertson effect from sun light fails by a factor of 40 to explain the 
observed rate; the effect from earthlight is even smaller. While a factor of 40 is close enough to give 
one pause, there are no parameters to adjust to try and nudge the result up to the required value. 

We rule these three mechanisms out. 

Resonance with the gravitational field and drag from interplanetary dust appear to be at least 
two orden-of-magnitude too small. However, these mechanisms contain “adjustable parameters" 

(the magnitudes of the gravity field coefficients and the density of the dust, for example) which 
might be increased to give the desired value. Also, we have used a simple linear theory in treating 
the resonance and ignored other complicating perturbations. Still, these two mechanisms appear to 
require unreasonable assumptions to attribute the force to either one of them. Szebehely ( 1 980), 
however, suggests the gravitational field is in fact responsible for the secular decrease; but our results 
on atmospheric drag strongly contraindicate this. 

Of the radiation effects the Schach effect may be ruled out. It operates only when part of the 
orbit is in shadow, while the secular decrease continues whether the orbit is in full sunlight or shadow. 
Further, it tends to increase the semimajor axis rather than decrease it. Terrestrial radiation pressure 
OZ the Yarkovsky effect do not seem to be the cause of the orbital decay. Terrestrial radiation 
pressure can increase and decrease the semimajor axis over a short period of time, while the Yarkov- 
sky effect should be absent in the early part of Fig. 1 when the orbit is in full sunlight. But their 
characteristic signature on the orbit is absent from Fig. 1 . However, the treatment of terrestrial 
radiation pressure and the Yarkovsky effect are the least satisfactory of all of those given here, the 
former because of its complexity and the latter because of its possibly large magnitude. More re- 
search is needed into their effect on satellite orbits. 
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Atmospheric drag seems to be the most likely cause for die secular decrease in the semimajor 
axis of Lageos’s orbit. Neutral hydrogen drag accounts for about 10 percent of the observed decrease. 
Charged particle drag from H+, He + , and 0 + accounts for at least another 61 percent. This makes a 
grand total of 71 percent of the observed decay which can be attributed to atmospheric drag. This 
is dm largest effect found so far with dm proper signature on dm semimajor axis. It is not hard to 
imagine that if this mechanism explains at least 71 percent, then it in fact accounts for 100 percent 
of the observed decrease. The argument is that our estimate of charged particle drag was somewhat 
conservative, so that in fact it contributes more than 61 percent of the total drag; and that our esti- 
mate of the neutral particle drag entirely omitted helium. A neutral helium component helps ex- 
plain the helium problem. 

In conclusion, it appears that atmospheric drag from a combination of charged and neutral par- 
ticles adequately explains the secular decrease of 1.1 mm day -1 in the semimajor axis of Lageos's 
orbit. 
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